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Abstract. Consider a N X n random matrix Y n = (Y?j) where the entries are given 
by Yg = zM^dM-Xg., the X§ being centered i.i.d. and a : [0,1] 2 -> (0, oo) being a 
continuous function called a variance profile. Consider now a deterministic N xn matrix 
A n = f-A-yJ whose non diagonal elements are zero. Denote by S„ the non-centered matrix 

Y n + A n and by A 1 " A n — min(7V, 71). Then under the assumption that lim^ — too — — c > 
and 

JVAn 



1=1 * ' 



where H is a probability measure, it is proven that the empirical distribution of the 
eigenvalues of S n Sj converges almost surely in distribution to a non random probability 
measure. This measure is characterized in terms of its Stieltjes transform, which is ob- 
tained with the help of an auxiliary system of equations. This kind of results is of interest 
in the field of wireless communication. 

Resume. Soit Y n = (Y") une matrice N X n dont les entrees sont donnees par Yfi = 
^C'/^j/") xVj , les A™, etant des variables aleatoires centrees, i.i.d. ct ou a : [0, l] 2 — » 
(0, oo) est une fonction continue qu'on appellera profil dc variance. Considcrons une ma- 
trice dctcrministc A n = f-A^J de dimensions N X n dont les elements non diagonaux 
sont nuls. Appelons T, n la matrice non centree definie par £ n = Y n + A n et notons 
N An = min(A, n). Sous les hypotheses que lim™—^ = c > et que 

JVAn 

I— 1 s ' 

ou H est une probabilite, on demontre que la mesure empirique des valeurs propres de 
E n Ejj converge presque surement vers une mesure de probabilite detcrministe. Cette 
mesure est caracterisee par sa transformed de Stieltjes, qui s'obtient a l'aide d'un systcmc 
d'equations auxiliaire. Ce type de resultats presente un interet dans le champ des com- 
munications numeriques sans fil. 



Key words and phrases: Random Matrix, empirical distribution of the eigenvalues, Stielt- 
jes transform. 
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1. Introduction 

Consider a N x n random matrix Y n = (Yy) where the entries are given by 

F „ = ^UM X , (L1) 
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where a : [0,1] x [0,1] — > (0, oo) is a continuous function called a variance profile and the 
random variables X™- are real, centered, independent and identically distributed (i.i.d.) with 
finite 4 + e moment. Consider a real deterministic N x n matrix A„ = (A^) whose non- 
diagonal elements are zero and consider the matrix E„ = Y n + A n . This model has two 
interesting features: The random variables are independent but not i.i.d. since the variance 
may vary and A n , the centering perturbation of Y n , though (pseudo) diagonal can be of full 
rank. The purpose of this article is to study the convergence of the empirical distribution 
of the eigenvalues of the Gram random matrix E n Ej (Ej being the transpose of E n ) when 
n — > +oo and N — ► +oo in such a way that — — - > c, < c < +oo. 

The asymptotics of the spectrum of N x N Gram random matrices Z n Z^ have been widely 
studied in the case where Z n is centered (see Marccnko and Pastur ^3] , Yin j2J , Silverstein 
et al. [HE], Girko |H1 E| , Khorunzhy et al. [Ej, Boutet de Monvel et al. gj, etc.). For 
an overview on asymptotic spectral properties of random matrices, see Bai |2]. The case of 
a Gram matrix Z n Z^ where Z n is non centered has comparatively received less attention. 
Let us mention Girko ([§], chapter 7) where a general study is carried out for the matrix 
Z n = (W n + A n ) where W n has a given variance profile and A n is deterministic. In [jjj, it is 
proved that the entries of the resolvent (Z n Z^~ zl)^ 1 have the same asymptotic behavior as 
the entries of a certain deterministic holomorphic N x N matrix valued function T n {z). This 
matrix- valued function is characterized by a non linear system of (n + N) coupled functional 
equations. Using different methods, Brent Dozier and Silverstein jS] study the eigenvalue 
asymptotics of the matrix (R n +X n )(R n +X n ) T in the case where the matrices X n and R n are 
independent random matrices, X n has i.i.d. entries and the empirical distribution of R n Rn 
converges to a non-random distribution. It is proved there that the eigenvalue distribution 
of (R n + X n )(R n + X n ) T converges almost surely towards a deterministic distribution whose 
Sticltjcs transform is uniquely defined by a certain functional equation. 

As in the model studied in this article, i.e. E„ = Y n + A„, is a particular case of 
the general case studied in ($.., chapter 7, equation K-j) for which there exists a limiting 
distribution for the empirical distribution of the eigenvalues. Since the centering term A„ is 
pseudo-diagonal, the proof of the convergence of the empirical distribution of the eigenvalues 
is based on a direct analysis of the diagonal terms of the resolvent (E„E^ — zl)^ 1 . This 
analysis leads in a natural way to the equations characterizing the Stieltjes transform of 
the limiting probability distribution of the eigenvalues. In the Wigner case with a variance 
profile (matrix Y n and the variance profile are symmetric) , let us mention the recent work of 
Anderson and Zeitouni [Q where the asymptotics of the spectum and central limit theorems 
are investigated by means of systematic combinatorial enumeration. 

Recently, many of these results have been applied to the field of Signal Processing and 
Communication Systems and some new ones have been developed for that purpose (Silver- 
stein and Combettes [T7], Tse et al. HUES], Debbah et al. 0, Li et al. [T3], etc.). The 
issue addressed in this paper is mainly motivated by the performance analysis of multiple- 
input multiple-output (MIMO) digital communication systems. In MIMO systems with n 
transmit antennas and N receive antennas, one can model the communication channel by a 
N xn matrix H n = (ii]J) where the entries if™- represent the complex gain between transmit 

antenna i and receive antenna j. The statistics C n = — logdet(/„ + H "f n ) (where H* is the 
hermitian adjoint and a 1 represents the variance of an additive noise corrupting the received 
signals) is a popular performance analysis index since it has been shown in information the- 
ory that C„ is the maximum number of bits per channel use and per antenna that can be 
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transmitted reliably in a MIMO system with channel matrix H n . Since 



C n =-5>g(l + ^), 



where (^k)i<k<N are the eigenvalues of H n H*, the empirical distribution of the eigenvalues 
of H n H* gives direct information on C n (see Tulino and Verdu j20] for an exhaustive review 
of recent results). For wireless systems, matrix H n is often modelled as a zero-mean Gaussian 
random matrix and several articles have recently been devoted to the study of the impact 
of the channel statistics (via the eigenvalues of H n H*) on the probability distribution of C„ 
(Chuah et al. [B], Goldsmith et al. |10 j . see also |20 | and the references therein). Of particular 
interest is also the channel matrix H n = Fn(Y u + A n )F^ where Fk = (Fp q )i< p , q <k is the 
Fourier matrix (i.e. F£ q = exp(2CT ^ p ~ 1 '^ q ~ 1 ^ )) and the matrix Y n is given by (|1 . 1|> (see 
|20j . p. 139 for more details). The matrices H n and E n having the same singular values, we 
will focus on the study of the empirical distribution of the singular values of £„. Moreover, we 
will focus on matrices with real entries since the complex case is a straightforward extension. 

In the sequel, we will study simultaneously quantities (Stieltjes kernels) related to the 
Stieltjes transforms of and £^£„. Even if the Stieltjes transforms of S„S^ and 
S^S„ are related in an obvious way, the corresponding Stieltjes kernels are not, as we shall 
see. We will prove that if N/n ► c > (since we study at the same time S„S^ and 



£jj£„, we assume without loss of generality that c < 1) and if there exists a probability 
measure H on [0, l]xK with compact support such that 



where T> stands for the convergence in distribution, then almost surely, the empirical distri- 
bution of the eigenvalues of the random matrix Yj n T^ (resp. E^E„) converges in distribution 
to a deterministic probability distribution P (resp. P). The probability distributions P and 
P are characterized in terms of their Stieltjes transform 





and 





l + C J CT 2 (- ,Cu)d7T z 



A 



H(du, dX) 



gdn z 




g(cu, A) 



H(du, dX) 



z(l + cj <T 2 (-,cu)dir z ) + 



A 



1+7 er 2 (u,-)d*z 




where the unknown parameters are the complex measures tt z and tt z . Both equalities stand 
for every continuous function g : H — > R, where Kc [0, 1] x 1 is the compact support of H . 
Then, this system admits a unique pair of solutions {jr z [dx, dX), n z (dx, dX)). In particular, 
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tt z is absolutely continuous with respect to H while tt z is not (see Remarks 12.51 12.61 and 12.71 
for more details). The Stieltjes transforms / and / are then given by 

f(z)= / Tr z (dx,dX) and f(z) = I Tr z (dx,dX). 

J[QS]xR J[0,l]xR 

The article is organized as follows. In Section |3 the notations and the assumptions are 
introduced and the main results (Theorem 12.31 and Theorem 12. 4|) are stated. Section [3] is 
devoted to the proof of Theorem 12.31 Section 0] is devoted to the proof of Theorem l2.4l In 
Section |3J we state as corollaries of Theorems 12 . 31 and 12 .41 two simpler results which already 
exist in the literature. We also give details for the statement of the results in the complex 
setting (Section |OJ). 



2. Convergence of the Stieltjes Transform 



2.1. Notations and Assumptions. Let N = N(n) be a sequence of integers such that 

nm = c. 

n— too TL 

Consider a N x n random matrix Y n where the entries are given by 
where XJ 1 - and a are defined below. 

Assumption A-l. The random variables (A™- ; l<i<N,l<j<n,n>l) are real, 
independent and identically distributed. They are centered with E(X™-) 2 = 1 and satisfy: 



3e>0, E\X"\ 4+e < oo. 



where E denotes the expectation. 

Assumption A-2. The real function a : [0, 1] x [0, 1] — > R is such that a 2 is continuous. 
Therefore there exist a non-negative constant cr ma x such that 

V(x,y) g [0,1] 2 , 0<<r a (x J »)<o^«<oo. (2.1) 

Remark 2.1. As long as a 2 remains continuous, a can vanish on portions of the domain 
[0,1] 2 . 



Denote by vax(Z) the variance of the random variable Z. Since var(Y^) = a 2 (i/N, j/n) / n, 
the function a will be called a variance profile. Denote by S zo ( dz) the dirac measure at point 
Zq. The indicator function of A is denoted by 1a(x). Denote by Cb(X) (resp. Cb(X;C)) 
the set of real (resp. complex) continuous and bounded functions over the topological set X 
and by ||/|| 

co — su-P^ga 1 |y(x)|, the supreimim norm. If ?C is compact, we simply write C(*-f) 
(resp. C(X;C))) instead of Cb(X) (resp. C b (X;C))). 

We will denote by —> the convergence in distribution for probability measures and by — ► 
the weak convergence for bounded complex measures. 

Consider a real deterministic N x n matrix A„ = (A^) whose non-diagonal entries are 
zero. We will often write Ay instead of Ay. We introduce the N x n matrix S n = Y n + A n . 
For every matrix A, we will denote by A T its transpose, by Tr(A) its trace (if A is square) 



ASYMPTOTIC DISTRIBUTION OF THE EIGENVALUES OF A GRAM MATRIX 



5 



and by F AA , the empirical distribution function of the eigenvalues of AA T . Since we will 
study at the same time the limiting spectrum of the matrices and S^E„, we can 

assume without loss of generality that c < 1. We also assume for simplicity that N < n. 

We assume that: 

Assumption A-3. There exists a probability measure H{ du, dX) over the set [0, 1] xK with 
compact support TL such that 

1 N 

i—l ^ ' 

Remark 2.2 (The complex case). Assumptions (A{Q, (A-|2J) and (A-|3J) must be slightly mod- 
ified in the complex setting. Related modifications are stated in Section 15.31 

When dealing with vectors, the norm || • || will denote the Euclidean norm. In the case of 
matrices, the norm || • || will refer to the spectral norm. 

Remark 2.3. Due to (A-|3J), we can assume without loss of generality that the A™-'s arc 
bounded for n large enough. In fact, suppose not, then by (A-|3J), jfJ2iLi^A 2 - ~* H\(dX) 
whose support is compact and, say, included in [0, if]. Then Portmanteau's theorem yields 
7? Ef ![o, K+S](^li) -> 1 thus 



#{A ii; Af^[0, K + S}} , 1 



= l-^E 1 M+^ A »)--0- (2-3) 



N N 

Denote by A„ = (A™ ) the matrix whose non-diagonal elements arc zero and set A"; = 
A"jl{(A".)2<K+a}. Then it is straightforward to check that ^(J- A 2 ^{du,dX) — ► 

H( du, dX). Moreover, if S n = Y n + A„ then 

||p ee- p ££-n W rank(S - S) W #{A ri , A| j[Q,K + 5]} (e) 

where (a) follows from Lemma 3.5 in |2J (see also ^5], Section 2), (b) follows from the fact 
that for a rectangular matrix A, rank(A) < the number of non zero entries of A and (c) 
follows from 1)2.3(1 . Therefore, F ss converges iff F ES converges. In this case they share 
the same limit. Remark 12.31 is proved. 

Remark 2.4. Due to Remark 12 .31 we will assume in the sequel that for all n, the support of 
jjJ23(±- a 2 ) 1S mcm dcd in a compact set JC C [0, 1] x M. 

Let C+ = {z e C, Im(z) > 0} and C v = {z e C+, |Re(z)| < Im(z)}. 

2.2. Stieltjes transforms and Stieltjes kernels. Let v be a bounded nonnegative mea- 
sure over M. Its Stieltjes transform / is defined by: 

JR A ~ z 

We list below the main properties of the Stieltjes transforms that will be needed in the 
sequel. 

Proposition 2.1. The following properties hold true: 
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(1) Let f be the Stieltjes transform of v, then 

- the function f is analytic over <C + , 

- the function f satisfies: \f{z)\ < jjSjj-j 

- if z G C+ then f(z) G C+, 

- if v(— oo,0) = then z G C + implies zf(z) G C + . 

(2) Conversely, let f be a function analytic over C + swc/i that f(z) G C + if z G C + 
and |/(^)| I Im(z) \ bounded on C + . Then, f is the Stieltjes transform of a bounded 
positive measure \i and //(K) is given by 

n(R) = lim -iyf(iy). 

y — '+00 

If moreover zf(z) G C + if z £ C + i/ien, /x(R _ ) = 0. 

(3) Let ¥ n and P be probability measures over R and denote by f n and f their Stieltjes 
transforms. Then 



faz G C+, f n {z) ► f(zj) 



Let ibe annxp matrix and let I n be the nxn identity. The resolvent of AA T is defined 

by 

Q(z) = (AA T -zln)- 1 = ( Pij (z)) x <. . j <„ , z G C\R. 
The following properties are straightforward. 
Proposition 2.2. Let Q(z) be the resolvent of AA T , then: 

(1) For all z G C+, ||Q(z)|| < (Im(z))- 1 . Similarly, \p i} {z)\ < (Im(z))- 1 . 

(2) The function h n (z) = ^Tr Q(z) is the Stieltjes transform of the empirical distribu- 
tion probability associated to the eigenvalues of AA T . Since these eigenvalues are 
nonnegative, zh n {z) G C + for z G C + . 

(3) Let I* be a n x 1 vector, then Im ^z^Q(z)^ r J G C+ for z G C + . 

Denote by .Mc('f) the set of complex measures over the topological set X. In the sequel, 
we will call Stieltjes kernel every application 

ft : C+ Mc(X) 

either denoted n(z,dx) or /j, z (dx) and satisfying: 



(1) Vg G Cb(A'), / gdfj, z is analytic 

(2) Vz G C+ V.g G C 6 (A-), 



< 



Im(z) 

(3) Vz G C+ V.g G Ct(A') and g > then Im (/ ff a> 2 ) > 0, 

(4) Vz G C+, V.g G C 6 (^) and g > then Im (z fgdfi z )> 0. 
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Let us introduce the following resolvents: 

Q n (z) = (S n E^ -zIn)' 1 = («t/( Z ))l<i,,j<iV> zeC+ > 
Q n (z) = (E^£„ - zln)' 1 = (qij {z)) x < id < n , zeC + , 
and the following empirical measures defined for z G C + 

1 N 

L n z {du,dX) = ^ S^^(du,d\), (2.4) 

i=l 
f " 

L"(du,dX) = - ^2qu(z) S/i_ A 2 j(du,d\) 

+ (- S±(du) (g> Sa(d\) J l {A r<„ } , (2.5) 

where ® denotes the product of measures. Since qu(z) (resp. qa{z)) is analytic over C + , 
satisfies |fe(z)l ^ (I m ( 2: )) _1 an( i mm (I m (<7ii(z))) I m { z 1u( z ))) > Oj £™ (resp. L n ) is a Stieltjes 
kernel. Recall that due to Remark |2. 41 L n and L n have supports included in the compact 
set JC. 

Remark 2.5 (on the limiting support of L n ). Consider a converging subsequence of L™, then 
its limiting support is necessarily included in TL. 

Remark 2.6 (on the limiting support of L n ). Denote by H c the image of the probability 
measure H under the application (u, A) i— > {cu, A), by TL C its support, by TZ the support of 
the measure l[ c l ](u) du ® 5o(dX). Let TL = TL C U 7?.. Notice that 7i is obviously compact. 
Consider a converging subsequence of L™, then its limiting support is necessarily included 
in TL. 

2.3. Convergence of the empirical measures L™ and L™. 

Theorem 2.3. Assume that H is a probability measure over the set [0, l]xK with compact 
support TL. Assume moreover that (A\gj) holds. Then the system of equations 

f 9dn * = I . r >, ,w |%v^ a H ( du ^) (2-6) 

J J -z(l + / a 2 (u, t)n(z, dt, dO) + 1+e/ ga(t|CM)7r(Z|dt|dC) 



gd7r z = c / — — - — t H(du,dX) 

-z(l + cj a 2 ft, CM )tt(z, dt, dfl) + i+ / ^( M ,t)*(z,dMC) 

+ (i-c)/' „ p#a , 



{1 + c Jo- 2 (t, u)n(z,dt,d()) 

(2.7) 

where \2. b]) and \2. 7| ) hold for every g g C(TL), admits a unique couple of solutions (tt(z, dt, dX), n(z, dt, dX)) 
among the set of Stieltjes kernels for which the support of measure -k z is included in TL and 
the support of measure n z is included in TL. 

Moreover the functions f(z) = J dix z and f(z) = J d~k z are the Stieltjes transforms of 
probability measures. 
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Remark 2.7 (on the absolute continuity of tt z and tt z ). Due to (|2.6|) . the complex measure 
tt z is absolutely continuous with respect to H. However, it is clear from l|2.7|) that tt z has 
an absolutely continuous part with respect to H c (recall that H c is the image of H under 
(u, A) i— > (cu, A)) and an absolutely continuous part with respect to 1[ Cj1 ](m) du Cg) <5o(^A) 
(which is in general singular with respect to H c ). Therefore, it is much more convenient to 
work with Stieltjes kernels tt and tt rather than with measure densities indexed by z. 

Theorem 2.4. Assume that (A^j), (A\^jj and (A\$) hold and denote by tt and tt the two 
Stieltjes kernels solutions of the coupled equations \2.b]) and \2. 7| ). Then 

(1) Almost surely, the Stieltjes kernel L n defined by \2-4\l converges weakly to it, that is: 

a.s. Vz e C+ L n z — tt z . 

n — >oo 

(2) Almost surely, the Stieltjes kernel L n defined by \2.5)) converges weakly to tt. 

Proofs of Theorems 12.31 and 12 . 41 are postponed to Sections and 

Corollary 2.5. Assume that (A^Tty, (A\^j and (A\$) hold and denote by tt and tt the two 

Stieltjes kernels solutions of the coupled equations \2. b}) and \2. 7| ). Then the empirical dis- 
tribution of the eigenvalues of the matrix S„S^ converges almost surely to a non-random 
probability measure P whose Stieltjes transform f(z) = J R+ is given by: 

f(z) = / TT z (dx,d\). 
Jh 

Similarly, the empirical distribution of the eigenvalues of the matrix converges almost 

surely to a non-random probability measure P whose Stieltjes transform f(z) is given by: 



f(z) = I TT z (dx,d\). 

J a 



Proof of Corollaru \2.5\ The Stieltjes transform of S n S^ is equal to j? J2i=i Qui 2 ) = / dL" 
By Theorem l2~4l 

a.s. Vz e C + , / dL n z ► / dn z . (2.8 



Since J dir z is the Stieltjes transform of a probability measure P by Theorem 12. 31 eq. (|2.8J) 
implies that F SnS » converges weakly to P. One can similarly prove that F s ™ s ™ converges 
weakly to a probability measure P. □ 
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3. Proof of Theorem 12.31 
Wc first introduce some notations. Denote by 
D(tt z ,tt z )(u, A) = -z ^1 + J a 2 (u,t)Tr(z,dt, d()^J + 

d(ir z )(u) = 1 + c / a 2 (t,cu)ir(z,dt,d(), 



1 + c J a 2 (t, cu)ir(z, dt, d() 



D(tt z , tt z )(u, A) = — z I 1 + c / cr 2 (t, cw)7r(z, di, d£) ) + 



1 + f a 2 (u,t)Tr{z,dt,d()' 



d(jt z )(u) = 1 + / cr 2 (u,t)Tr(z,dt, d(), 
k(7T 2 )(m) = — z 1 + c / <r 2 (t 7 u)ir(z,dt 7 d£) 



Let v be a complex measure over the set H (recall that H is compact by (A-|5J|) then we 
denote by ||^||tv the total variation norm of v, that is 



IIHItv = W\(H) 



sup 



fdv 



, /eC(H;C), ll/Hoo < 1 



3.1. Proof of the unicity of the solutions. Notice that the system of equations 
and 12.7fl remains true for every g £ C(H;C) (consider g = h + ik) and assume that both 
("zjTTz) and (p z ,Pz) are pairs of solutions of the given system. Let g € C(Tt;C), then (|2.6|l 
yields: 

/" [ zg(u,X) f cr 2 (u,t)(n(z 1 dt,d() - p{z,dt,d()) 

gdir z - gdp z = / ^— - — H(du,dX) 

J J D{tt Zi tt z ) x D{p Zl p z ) 

cX g(u, X) J a 2 (t, cu)(p(z, dt, d() - ir(z,dt,d()) 

; r ; — — r tl (CLU, UA I 

D(w z ,tt z ) x D(p Zl p z ) x d(ir z ) x d(p z ) 
and 

dH 



g dir z — g dp 



< kkmaxllffllooll^ ~Pz\\u 



CCTmax llfllloo \Wz ~ Pz\\tv 



\D(tt z ,tt z ) X D(p z ,Pz)\ 

XH{du, dX) 



\D(tt Zi tt z ) x D(p z ,p z ) x d(7r z ) x d{p z )\ 



If one takes the suprcmum over the functions g G C(H; C), ||flf||oo < lj one gets : 
IKz - /3z||tv < "Ikz - Pzjltv + /3||7r z - /5 2 ||tv 

where 

. _ 2 f XH(du,dX) 
a = a(TT,p,ir,p) = ccr max / — — — — — — — — -r 

7 \D{it z ,h z ) x D{p z ,p z ) x d(7r z ) x d(p z )| 

\D(tt Zi it z ) x D(p a ,p z )| 
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Similarly, lt!T7l) yields: 
g dir z - g dp z 



/ czg(u,X) J a 2 it, cu){w(z,dt, d£) - p(z, dt, dp) g 



-D(7r z ,7r z ) x D(p z ,p z ) 
Xg(u, A) / cr 2 (ii, i)(7r(z, di, d() — p( z > dt, d()) 

D(tT z ,TT z ) X D{p z ,p z ) X J(7T Z ) x d(p 2 ) 



H(du,dX) 



_ : l c z g(u, 0) / a 2 (t, u){tt(z, dt, d() - p(z, dt, dp) ^ 
k(tt z ) x n(p z ) 



and 



gdir z - / g dp 



< c2cr maxkl ll9l|oo|k 2 -p 2 || tv 
T" C (7 max 

II.9II ooH"^^ Az||tv 



dH 



b{-K z ,Tt z ) X D(p z ,p z ) 

XH(du, dX) 



b(n z ,Tr z ) x D(p z ,p z ) x d(?r z ) x d(p z ) 



+ (1 - c)ccr^ ax |^| ||g||oo||7r 2 - p 2 || tv 



r/(( 



|k(7T z ) X K(pz)\ 

As previously, by taking the suprcmum over g G C(7i; C), ||g||oo < 1 ; wc get: 
||7Tz - (Ozlltv < "Iks - PzWtv + PW^z - PzWtv 



where 

a = ct(n, p, 7r, p) 



C<J mav z c 



dH 



= /3(w, p, 7f , p) = cal 



D(it z ,tt z ) x D(p z ,p z ) 
XH(du,dX) 



(1-c) 



du 



\k(tt z ) x n(p z )\ 



£>(tt z ,tt z ) x D(p z ,p z ) x d(7r 2 ) x d{p z ) 



We end up with the following inequations: 



Kz - p z || tv < ol\\tt z - p z \\ tv +/3\\tt z - p z \\ tv 

\TTz - Pz\\tv < <5|kz - Pz\\tv + P\\Kz - Pz\\tv 



(3.1) 



Let us prove now that for z G C v with lm(z) large enough, then a < i. 

Since 7r and fr are assumed to be Stieltjes kernels, Im (z J (T 2 (u, i)7r(z, dt, dQ) > and 
Im (J a 2 (t,cu)n(z,dt,dC)) > 0. Therefore, Im(£>(7r 2 , 7r 2 )) < — Im(z) and hence |Im(£>(7r 2 , 7r z ))| > 
Im(z). Similarly, \Im(D(p z , p z ))\ > Im(z). Thus, 



< 



\D(Ttz,n z )xD{p z ,Pz)\ _ Im 2 (z)' 
Now consider zd(ir z ). As previously, Im(zd(7r 2 )) > Im(z). As |zd(7r 2 )| > |Im(zd(7r 2 ))|, this 



implies that 



< 



and 



\zd(TT z )\ — Im(z) 

argument holds for d(p z ) thus we get 



Wkr\ < Since * e C v , J£L_ < V2. The same 



2ccr 2 lax f XH{du,dX) 1 / \ i 

a < 2 < — i° r fm(^) large enough. 

Im (z) 2 
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With similar arguments, one can prove that 

/3<^f, a<f4r' /3<-^tfr f XH(dudX). (3.2) 
Im(z) lm(z) Im (z) J 

Therefore max(a, (3, a, /3) < 8 < \ for z G C v and Im(z) large enough where 9 does not 
depend on (71", fr, p, p). Therefore, the system (|3.1|1 yields 

— Pz\\tv — W^z — Pz\\tv = for z G C v and Im(z) large enough. 

Now take z G C + and g G C{TL). Since f gdir z and J gdp z (resp. J gd7r z and J gdp z ) are 
analytic over C + and are equal in C v for Im(z) large enough, they are equal everywhere. 
Since this is true for all g G C(H), tt z and p z (resp. tt z and p z ) are identical on C + . This 
proves the unicity. 

3.2. Proof of the existence of solutions. Let us now prove the existence of solutions to 
(|3.1|) . Define by recursion 

7r°(z, du, dX) = tt°(z, du, d\) = H(du, dX), 

z 

and 

' Jg(u 1 X)7rP(z,du,dX) = J ^0ZL^ H(du,d\) 
J g(u, \)n*(z, du, dX) =cj ^*>,_ L) H(du, dX) + (1 - c) £ d« 

for all g G C(H). It is straightforward to check that tt z (resp. tt z ) is a Stieltjcs kernel. 
Moreover, this remains true for n p and tt p by induction over p. As for the unicity, we can 
establish 

kf-^ritv <aikr 1 -^ _2 iitv+i9||7f5- 1 -7fr a iitv 
\n - Tff-^itv < fiikr 1 - ^r 2 ikv + ^ikr 1 - ^- 2 i! tv [ ' 

where a,a,f3 and (3 depend on (tt^ -1 , 7r^ _1 , t^ z ~ 2 , ^f" 2 )- As in 1)3. 2|) . one can prove that 

max(a, a, (3, 0) <0<— for z G C v and Im(z) large enough 

where 8 does not depend on (^p -1 , 7ff _1 , 7rf -2 , 7ff -2 ). Therefore, (7r§) and (fr£) are Cauchy 
sequences with respect to the norm || ■ || tv whenever z G C v and Im(z) is large enough. This 
yields the existence and unicity of kernels ir z and 7r 2 such that 

{ h^l > JgdTTz 

Vg G C(H), I r , 

whenever z G C v and Im(z) is large enough. 

Let g G C(H) be fixed. Recall that z t—> J gdir? is analytic on C + and that Vz G C + , 

\fgdnP\ < . Therefore, (/ gdTT^) p >o is a normal family. From every subsequence of 

(/ gdnP)p one can thus extract a converging subsequence (j g dit z I )M where M = M(p) 
such that 



VK C C + , isT compact sup 

z<£K 



J gdnf -Y{g){z) 



M—>oo 



0. 



where F(<7)(z) is analytic over C + . If z G C v and Im(z) is large enough, we know that 
J g dnP — > J gdir z . Therefore, T(g)(z) — J gdn z for z G C v and Im(z) large enough. 
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From this we can conclude that for all z G C + , every subsequence has the same limit, 
say T(g)(z) thus: 

Vz G C+, f gdn* ► r(.g)(z), (3.4) 

J p-+oo 

where T(f)(z) is analytic. Moreover, it is straightforward to prove that 

(1) T(ag + bh)(z) = aT(g) + bT(h), 

(2) \T(g)(z)\ < Mfy, 

(3) ha(V(g)(z)) > if g > and z G C+, 

(4) ha(zT(g)(z)) > if g > and z G C+. 

As W is compact and since the application g i— > r(<?)(z) defined for g £ C(W) is lin- 
ear (property ^) and continuous (property |2J) , the Riesz representation theorem yields the 
existence of a measure ■K z {du 1 dX) = n(z, du, dX) such that 



T(g)(z) ~ j g(u,X)ir(z,du,d\). 

Similarly one can prove that 

J gd%l^V{g){z) = Jgdn z . (3.5) 

Let us now prove that it and tt satisfy 1)2. 6fl and (|2.7[) . We first check that 1 

toeC + ,Vue [0,1], rf(7r z )(it) 7^ 0. (3.6) 

Indeed assume that for a given u, there exists zq G C + such that d(ir Zo )(u) = and 
consider the function $(z) = / a 2 (t,cu)ir(z,dt,d(). As d(7r z )(u) = 1 + $(2), we have 
Im($(zo)) = 0. Since Im(<I>(z)) is harmonic and non-negative over C + , the mean value 

property implies that Im(<&(z)) = over C + . By the Cauchy-Ricmann equations, Re(<i>(z)) 

2 

is therefore constant. But since |$(z)| < , we have $(z) = 0. This yields in particular 
d(jr Zo )(u) = 1 which contradicts d(ir Zo )(u) = 0. 

Due to (El, E3 and E3, one has D{ ^Jv ){uA) ^ P D{i J r){u , x y Since g(M)KA) ^ 
j— j^j, the dominated convergence theorem yields: 

r^t-p P w — 7-rH(du,d\) ► / — — -H(du,dX). 

D(lTz, 7Tz){u, A) J D(tt z , tt z )(u, X) 

On the other hand J g dit v z — > p J g dir z and 12. 6|) is established. One can establish Eq. (|2.7|l 
similarly. 

It remains to prove that f(z) — J dir z is the Stieltjes transform of a probability measure 
(one will prove similarly the corresponding result for f(z) = J dw z ). Recall that 

hn(f(z)) = Im ( ( diiA > [ ^H(du,dX)>0 

UV " \J J ~ J \D(n z ,n z )(u,X)\ 2 1 ' 

by (|2.6|) . Moreover, since < j^r^ , f(z) is the Stieltjes transform of a subprobability 

measure. It remains to check that lim y _> +00 iyfiiy) = — 1. Since 

T 2 



<7 (u,t)ir(iy,dt,d() 



9 

< and 

y 



a (t,cu)ir(iy,dt,dQ 



4l existe peut-etre un argument plus naturel base sur |d(7rf)| > A voir. 



ASYMPTOTIC DISTRIBUTION OF THE EIGENVALUES OF A GRAM MATRIX 



13 



and 

iy H{du, dX) 



iyf{iy) 



-iy(l + J a 2 (u, t)n(iy, dt, dQ) + 1+c/ g 2 (t|Cu)ff(<y|df|dC) 
the Dominated convergence theorem yields the desired result. Theorem 12.31 is proved. 

4. Proof of Theorem 12.41 
We first give an outline of the proof. The proof is carried out following three steps: 

(1) We first prove that for each subsequence M(n) of n there exists a subsequence 
Msub = M su b(n) such that for all z £ C + , 

£ Af.» b _!!!_ ^ and £M sub ^^- zj (41) 

n — >oo n — >oo 

where [i z and Ji z are complex measures, a priori random, with support included in 
H (Section EJ). 

(2) We then prove that z i— > \i z and z i— ► p, z are Stieltjes kernels (Section |4.2(l . 

(3) We finally prove that for a countable collection of z £ C + with a limit point, say 

C = {zpjpeN U {zoo} with Zp > x-oo ; 

the measures ijl z and /iz (which are a priori random) satisfy equations 1)2.6(1 and 1(2.7(1 
almost surely for all z £ C. Since C has a limit point in C + , analyticity arguments 
will yield: 

a.s. Vz £ C + , n z = ir z and £l z = tt z . 
Otherwise stated, 

a.s. Vz £ C + , L" ^> 7r 2 and L™ ^> n z 
which yields the desired result (Section (4. 3f l. 



4.1. Step 1: convergence of subsequences L* /sub and L* /sub . Let zq £ C + and let B = 

{z £ C, \z — zo\ < 8} C C + . Due to assumption (A-OJ) and to the fact that |<jzi(z)| < Im _ (z), 
Helly's theorem implies that for each subsequence of n there exists a subsequence M = M(n) 
and a complex measure fi Zo such that 

n — >oo 

Since L n is random, \i Zq is a priori random too but due to (A-I2J), its support is included 
in H. Let (zk,k > 1) be a sequence of complex numbers dense in C + , then by Cantor 
diagonalization argument, one can extract a subsequence from M, say M su b, such that 

Vfc£N, L%«* p Sk and L^f sub — - — ► p, Zk , 
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where [i Zh and fl Zk are complex measures, a priori random. Let g € CQC) and let z S 
There exists Zk such that I z — Zk I < e and 



(a) 



, dL M Bub („) _ / gdL M aab (m) 



< 



gdLf'^ n) 



g dL 



(f>) (c) 

Let n and m be large enough. Since £^ sub converges, (b) goes to zero. Since quiz) is analytic 
and since < Im _1 (;z), there exists K > 0, such that 

Mi > 1, Vz, z' close enough, \qu(z) — qu(z')\ < K\z — z'\. 

Thus max{(a), (c)} < A'HgHoolz — z^\. Therefore, (J gdL^ 1 ^) is a Cauchy sequence and 
converges to Q(g)(z). Since g i— > 0(.g)(;z) is linear and since |0(g)(z)| < Im _1 (z)||g|| 00 , Riesz 
representation's theorem yields the existence of /x 2 such that 0(<?) (z) = J gd/j, z (recall that 
the support of /i z is included in Ti which is compact). The convergence of L^ 1 "^ can be 
proved similarly and l|4.1|) is satisfied. The first step is proved. 



4.2. Step 2: the kernels \x z and fi z are Stieltjes kernels. Let us now prove that 
z i— ► J gd\i z is analytic over C + . Since |f gdL^ Isub \ < Im~ 1 (-2 ) 1 1 5 1 1 oo 5 from each subsequence 
of (J g dL^ sub ) , one can extract a subsequence that converges to an analytic function. Since 
this limit is equal to J gdfi z , the analyticity of z > J gd\i z over C + is proved. Since 
properties J3J and |J2J defining the Stieltjes kernels are satisfied by L™, the kernel \i z inherits 
them. Therefore, \i z is a Stieltjes kernel. Similarly, one can prove that \x z is a Stieltjes kernel. 
The second step is proved. 



4.3. Step 3: the kernels fi z and pt z are almost surely equal to w z and tt z . We will 
now prove that almost surely \x z and p, z satisfy equations (|2.6|) and (|2.7(l . 

In the sequel we will drop the subscript n from the notations relative to matrices, and 
the superscript n from A,™. Let e*i = (<5i(fc))i<fc< n and fi = (6i(k))i<k<N ■ For the sake of 
simplicity, S T will be denoted S. Consider the following notations: 



Matrix 


Y 


A 


s 






A' J ' 




77T 
~M 


zth row 


Vi- 




& 






A«/j 




Vi 


Matrix when ith row 
is deleted 








^(i,t) 











In particular, £j. = + A^e^ and = y.,- + Ajj/j for 1 < i < N . We will denote by Di 
and Aj the respectively n x n and N x N diagonal matrices defined by 



A = diag 



diag 




Finally, for 1 < i < N, we denote by -D(M) an( l ^(M) ^ ne ma trices that remain after deleting 
row i and column i from Di and Aj respectively. 



We can state our first lemma: 
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Lemma 4.1. Let z G C + be fixed. 

(1) The ith diagonal element qa(z) of the matrix (SS T — z/jv) -1 can be written : 

qU ^ Z) = z z T n jM' k )5.Jz) I A '- 1 A I c (2) I c (3) 

(4.2) 

where 1 < i < N and 
£ g = ( E « E « - € + fl>? (Ef^Ep) - zl) 1 



e (4) 



( E (v) E M C-Trf^v) ( E (M) E CM) - z/ ) ) 



£ S = ^( A (M) ( E M S CM)-* J ) J -1^(A?(ES T -z/)" 
Moreover almost surely, 



1 ^ 

Vfc,l<fc<5, lim -$>2 =0. (4.3) 



i=l 



(2) If 1 < £ < N then the ith diagonal element qu(z) of the matrix (S T S — zl n ) 1 = 
(SS T — zln)^ 1 can be written: 



i+iEI=i-' 2 (^,i)i«W+^+^ 
If N + 1 < i < n, then qu can be written: 



(4.4) 



w/iere 



efl = zTrtAUErZ-ziy^-zTrfAllE^-ziy^ 
if.i (s^Efij - zlj if, - Tr (d 2 (i i) (e^E^ - zlj 



£ -(4) 
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Moreover, almost surely 
for k - 



1 N 

1,4,5 lim -V ef 



1 " I 

for k = 2,3 lim - V gV „ = 0. 



(4.6) 



Proof of Lemma \4-1\ Since qu(z) = (SS T — zl)^ , this element is the inverse of the Schur 
complement of ^£(j)£j^ — zl^j in (SS T — z/) (see for instance QT], Appendix A). In other 
words 

?«(*) = (||fi.|| 2 -«-g.E5 ) (E (i) E5 ) -zJ)- 1 E w ^')" 1 . 
Using the identity 

we have 

1 



<&(z) = 



-z — zt 



- f Efe=i ^ 2 (&> I) - ^(^(i) - zl^ej + A ii£ g + eg + eg ' 

Similarly, we have 

- zir'ef = (Ef i} S (i) - zl)^ 



-z - zf) i .(£(ii)'£T. ) - zl) l fj. 

1 



ST 



(4) , ,(5) 



(4.7) 



And i|4.2[l is established. It is important to already note that since ffi. is the ith row of , 
Effi. = (while Eyj. = (0, • • • ,A it , ■ ■ ■ ,0)). If i < N, (|4.4jl can be established in the same 
way. If i > N + 1, then £.j is centered: There arc no more Aji and all the terms involving 
Ajj disappear in 1)4.4(1 . which yields 1)4. 5 [I . 



We now prove that 



-T e (1) 



N 



(4.8) 



One will prove similarly that ^^1=1 l%nl — * a - s - Denote by R n = (E^E(j) — zl) 1 = 
(pij). Since i?„ is symetric, ej„ = —z2iji.R n e[ and 



\yi.R n ef | 4 = I ^ l^fcPfeil 4 = X] YikiYikiY^YuzPta 



i Pk 2 i Phi Phi- 
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Denote by E^ n the expectation conditionnally to the er-algebra generated by R n . Since y^. 
and R n are independent and since EY^ — 0, we get: 

E^J^efl 4 = 2E Rn Y ?M 2 YZ\pu\ 2 

+E Rn ]T YlPk^uPu+^Y,^^ 

k,l;kjtl k 

< 4E(X^) 4 ^^|p fci | 2 | Pw | 2 =4E(Xg.) 4 %^ (E 
but E fc |Pfc t | 2 = ||i?„e 2 || 2 < ||i?„|| 2 < j^y. Therefore, 



Pki\ 



E 



.(i) 



4 |2 2 | 4 4E(Xg)V: 



/ max 



n 2 Im 4 (z) 



(4.9) 



Finally. 



1 N 

-Y 



><5> 



6 4 iV 
1 



< 



5 4 N 



-AY 



_(1)|4 



/ 1 TCP I (1) |4 ( b ) 1 

o kkjv n 



where (a) follows from Minkowski's inequality and (b) from (|4.9|l and Borel-Cantelli's lemma 
yields Eq. (|4~%|) . 



Let us now prove that 



N , 

-T M 2) 



(4.10) 



One will prove similarly that ^ X] l^inl; a? S \ £ i 4 l\ an d ^ S l^i^l S° to zero a - s - Denote by 
Xi. = (Xn, ■ ■ ■ ,Xi„) and write y*i. = Xi.Di. In particular. 

ft(£5j£ (i) - z/)" 1 ^ = £,A(S^S W - ziy'Djxl 

where 2J£. and Dj(E£sE(j) — zI)~ 1 Df are independent. Lemma 2.7 in [3] states that 



El^.Cff -TrC| p 



p ((E(X«)' 



Tr CC 



T\P/ 2 



+ E(X 4l ) 2p Tr(CC T ) p/2 ) (4.11) 



for all p > 2. Take p = 2 + e/2 where e is given by (A-fTJ and let C = D^E^E^ - zI)- 1 Df. 



Then 



i'7 



Vg > 1, Tr{CC T ) g < 



Therefore, lj4"TT) and i jlTT^I yield 



Im 2 «(; 



(4.12) 



where the constants K, K\ and K2 depend on the moments of Xn, on cr max and on Im(z). 
Thus 



E 



.(2) 



p < K\z\ p 



,1+6/4 ■ 



(4.13) 
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Finally. 



N ^ 



.(2) 



(<*) i 
< 



SpNp 



' N 

E( e i^i p )' 



\i=l 



< ± sup E|egr S * 
o p Ki<N n i+£ / 4 



where (a) follows from Minkowski's inequality and (b) from (|4.13|) . and Borel-Cantclli's 
lemma yields (|4. 101) . 
We now prove that 



1 N 

-Y 



.(3) 



2 = 1 



0. 



(4.14) 



One will prove similarly that ^ l^i ™l S oes to zero. Since E T E = S^E^j + Lemma 
2.6 in [ini yields: 



Tr 



((E T E- 2 7)- 1 -(E^ ) E w -. ! 7)- 1 



A- 



< 



nlm(z) ' 



In particular, 



.(3) 



< 



nlm(z) 



(4.15) 



which immediatly yields (|4. 14|l . 
We finally prove that 



i N , 

-T M 5) 



i=i 



(4.16) 



One will prove similarly that ^ goes to zero. Write 



.(5) 



e i,n = Tr Afi t {) 



S fi,i) - Zl 



Tr A? (EE T -zI 



As for e^, one can prove that 



TrAf,,, E«. fl E 



2/ 



v (i,i) ^(i,») Zj (t,i) 

by applying Lemma 2.6 in |16| . Let 



TrA (M) ( E (i) S (i)-* J 



< 



nlm(z) 



= Tr A (M) (E (i) E (i) - 



Tr A, (EE - zJ 
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By applying to SS T — zl the identities relative to the inverse of a partitioned matrix (see 
[TT] . Appendix A), we obtain: Tr A^ (EE T - zl)~ l = *i + * 2 + *3 where 



#2 = 



i£r 



#3 = -- 



- 2 (lv4) 



^-^-4.(E^S (i) -z/)-^ ' 
In particular, Kj ; „ = — \E r 2 — ^3- We have 



|* 2 | = 



^.(Ef E w -z/)-^ 



< 



KM) I 



Let S(j) = X^i* v i u i v J be a singular value decomposition of where i/;, ui, and are 
respectively the singular values, left singular vectors, and right singular vectors of E(j). Then 



(e«e^-^)" s (i) eT =E 



and 



Im (z + - z/)- 1 ^) - Im(z) 1 + £ 



1=1 Wf- Z \ 



( iv-i^Lr^ 2 



i=i i 



2 I 2 



V 

As a consequence, 1*21 < "if^i^j)- Furthermore, since Im fz|i.(£^£(j) — > 

by Proposition (3), we have |* 3 | < ^f^i^y- Thus, ef^ < ^^j^y, which immedi- 
atly yields (|4.16|) . Lemma f4. II is proved. □ 



Recall notation D introduced at the beginning of Section 3: 
D(tt z , ir z )(u, A) =-z ( 1+ / <r 2 (u,t)n(z,dt,d() ) + 



1 + c J a 2 (t, cu)ir(z, dt, d() ' 



Corollary 4.2. Let z G C + 6e fixed. Then almost surely 

N , A 



V.g e C(/C), lim 



g(A%,i/N) 



= 0. 



(4.17) 
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Proof of Corollary \4-<H Following the notations D and d introduced at the beginning of 
Section [31 we introduce their empirical counterparts: 



D n (u) 



k=l 

Z 



£ (4) +e (5) 

i.n 1 i.n 



^2<y 2 (u, - J + 

fe=l ^ n ^ 



+A-e (1) +e (2) +e (3) 



A 2 

ft 

d"(u) 



Since q u = {D n {i/N))- 1 by 63 and (S^S (i) - z/)^ 1 = (i/JV)) -1 by gTJ), Propo- 

sition El (1) yields: 



1 1 

< — and 



i < M 



|£>"(i/JV)| ~ Im(z) ' |d"(i/iV)| ~ lm(z) 
On the other hand, since /i z and /2 Z are Stieltjes kernels, we have: 



1 



< 



1 



and 



< 



Therefore, 



|D(/2„/i,)(i/iV,A^)| " Im(z) —* \d(fi z )(i/N)\ " Im(^ 

1 -z(j<?(i/N r )dL«-fo*(i/N,-)dii z ) 



(4.18) 



(4.19) 



D(fi„ A*) 



D-*(i)xD(fc )fe )(iA|) 

A| (| / a 2 (., »/n) dLg - c/ a 2 (-, »/n) 
(£) x d n (£) x Z> (±) x (i, A?.) 

Al^g+jg) 

+ dO*,) (£) x d™ (i) x D» (i) x D(£ z ,/x.) (i7, A 2 ) 

Recall that the A^'s are assumed to be bounded (say |A^| < K). Due to (|4.18|) and 
Jill , we get: 



D(/2 Z ,^)(^,A|) 
(i) 



< 



W(z) 



a 2 (i/N,-)dL™~ a 2 (i/N,-)dfl 



I(i,n) 



K 



+ 



.(2) 



.(3) 



|z| 2 if 2 
Im 4 (z) 



Im 2 (z) 
/ 

N 



a 2 (-,i/n)dL" z - c / a 2 {-,i/n)dn 



.(4) 



.(5) 
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In order to prove sup i<N n) -> 0, recall that C([0, l] 2 ) = C([0, 1]) ® C([0, 1]). In 
particular, Ve > 0, there exists k G N, gz G C([0, 1]) and ft.; G C([0, 1]) for I < k such that 



su Pl 



(J 2 (x,t)-Y:^i9l(x)hi(t) 



< e. Therefore, 



cr 2 (i/N,-)dL r z 



a 2 {i/N,-)djl z 



< sup 



a 2 (x,-)dL* 



a 2 (x,-)dp, z 







which implies that sup i<A r goes to zero. One can prove similarly that sup i<Ar J(i, n) 

goes to zero. Therefore, 



N N 



<?(A?^/iV) 



A^ 



i=i 



< |z|l L g|l °° B upf(«,n)+ NU 



Im (z) i<N 



InT(z) 



N N N 

i=l 



1=1 



-(3) I 



i=l 



Izp^Hslloo / 1 ^ (4), , 1 V^i (5), 



and l|4.17[) is proved with the help of Lemma 14.11 



□ 



We now come back to the proof of the third step of Theorem 12.41 For simplicity, we will 
denote by n* = M su b(n) where M su b is defined previously, by A^* = N(n*). 



3(A") 



A direct application of the Dominated convergence theorem yields that (A, it) . - - ^ u ^ 
is bounded and continuous therefore (A-|3J) yields 

^£ ^ (A ^:>^ — f „^± „ w («»> 



Moreover, 



^X>(A„M*)(*/JV*,A£) ™ 7 D(A»,M»)(«,A) 



i N ' r 

l^I^A*. *'/#*)«« ► hdfi, 



(4.21) 



Consider now a countable set C with a limit point. Since C is countable, l|4.17[) holds almost 
surely for every z G C and for every g G C{JC). Thus l|4.20|l and H4.21f> yield that fj, z and jl z 
satisfy l|2.6|) (and similarly l|2.7[l ) almost surely for all z G C. 

Since \x z and fl z are Stieltjes kernels, one can easily prove that z i— > J jjr~^ ^ \ dH is 

analytic over C + . Therefore, by l|2.6[) . the two analytic functions z i— > Jgdfx z and z i— > 



rdH coincide almost surely over C which contains a limit point. They must be 



equal almost surely over C + . Therefore fi z and £l z satisfy l|2.6[) (and similarly H2.7|) l almost 
surely for all z G C + . 

Since \i and jl are Stieltjes kernels satisfying almost surely l|2.6[) and (|2.7|l . they must be 
almost surely equal to the unique pair of solutions (tt, fr) by Theorem 12.31 In particular, /x 
and /i are almost surely non-random. Thus for every subsequence M = M(n), 



a.s. 



Vz G 



M w 



and L 



M w 



Therefore, the convergence remains true for the whole sequences L™ and L™. Theorem 12.41 
is proved. 
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5. Further Results and Remarks 

In this section, we present two corollaries of Theorems 12.31 and 12.41 We will discuss the 
case where A n = and the case where the variance profile cr(x, y) is constant. These results 
are already well-know ( g| E| ) • 

5.1. The Centered case. 

Corollary 5.1. Assume that (04-QJ) and (A\^ hold. Then the empirical distribution of the 
eigenvalues of the matrix Y n converges a.s. to a non-random probability measure P whose 
Stieltjes transform f is given by 



f{z) = / -Kz(dx), 

J [0,1] 

where tt z is the unigue Stieltjes kernel with support included in [0, 1] and satisfying 



s(«) 



V ff eC([0,l]), gdn z = -j -du. (5.1) 

Remark 5.1. In this case, one can prove that tt z is absolutely continuous with respect to du, 
i.e. TT z (du) = k(z,u)du where z <— > k(z,u) is analytic and u i— > k{z,u) is continuous. Eq. 
I|5.1(l becomes 

Vug [0,1], VzeC+, k(u,z) = . (5.2) 

Eq. (|5.2(l appears (up to notational differences) in [H] and in 0] in the setting of Gram 
matrices based on Gaussian fields. 

Proof. Assumption (A-|3J) is satisfied with A^ = and H(du,dX) = du (g) <5o(A) where d« 
denotes Lebesgue measure on [0,1]. Therefore Theorems 12.31 and 12.41 yield the existence of 
kernels ir z and n z satisfying (|2.6|l and i(2.7() . It is straightforward to check that in this case 
tt z and 7r 2 do not depend on variable A. Therefore 1)2.6(1 and 1)2.7(1 become: 

gdir z = [ g 2 ^, w ,,^ du ( 5 - 3 ) 

J — z(l + J o- z (u,t)Tr{z, at)) 



and 



g[cu) du 



(i, cu)7t(z, dt)) 



+ (1 ~ C) /[ci] -z(l + cja^t,u)n(z,dt)) 

3>1] -z(l + c/ a 2 (t,u)^(z,rft)) dU ' (5 ' 4) 

where g S C([0, 1]). Replacing / a 2 (u,t)ir(z,dt) in 1(5.3(1 by the expression given by 1(5.4(1 . 
one gets the following equation satisfied by ir z (du) : 

g( u ) 

+ f o 2 (u,t) ,r 

^ J l+c [ cr 2 (s.t)n(z,ds) U ' L 



gdn z = / 57 — -t du 

□ 
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5.2. The non-centered case with i.i.d. entries. 

Corollary 5.2. Assume that f-A-yJ) and (A\$) hold where o~(x, y) — a is a constant function. 
Assume moreover that jf^2iLi°~A 2 - ~~ * H\(dX) weakly, where H\ has a compact support. 
Then the empirical distribution of the eigenvalues of the matrix Yi n T^ converges a.s. to a 
non-random probability measure P whose Stieltjes transform is given by 

f(z) = f M^l (5 5) 

J - 2 (l + «7 2 /(z)) + (l-c)^ + TT ^ 7R - 

Remark 5.2. Eq. (|5.5|l appears in (Kj in the case where S„ = aZ n + R n where Z n and R n are 
assumed to be independent, Zfj = — jS-, the a,j being i.i.d. and the empirical distribution 

of the eigenvalues of R n R„ converging to a given probability distribution. Since R n is not 
assumed to be diagonal in the results in jS] do not follow from Corollary 15. 21 

Proof. One can build a sequence (i/n, Af 4 ) such that ^ %/n.A 2 ) ~* du ® H\(dX). 
Therefore (A-|3) is satisfied with H(du, dX) = du® H^(dX) and Theorems 12.31 and 12 . 41 yield 
the existence of kernels ir z and tt z satisfying l|2.6[) and 12.7|l . It is straightforward to check 
that in this case ir z and tt z do not depend on variable u. Equation (|2.6|l becomes 

gdir z = [ - x H{du,dX) 

Let g(u, A) = 1, then i|2.6[) becomes 

/(*) = / — 2^7~\\~i x H ^ ( 5 ' 6 ) 

J -z(l + a 2 f(z)) + 1+ J f(z) 

Denote by f n (z) = ^Y,i Qa( z ) and b y = ^ <fe( z ) = tT Tr ( S n S » _ zl)' 1 . Since 

f n (z) = iTr(£„E£ - zl)- 1 and f n (z) = iTr(S^E„ - zl)' 1 (recall that N < n), we have 
fn(z) = £/„(*) + (1 - f ) . This yields /(z) = c/(z) - i=2. Replacing />) in 

by this expression, we get l|5.5|l . □ 

5.3. Statement of the results in the complex case. In the complex setting, Assump- 
tions (A-^)-(A-|2J| must be modified in the following way: 

Assumption A-l. The random variables (X^ ; l<i<N,l<j<n,n>l) are complex, 
independent and identically distributed. They are centered with E|A™ | 2 = 1 and satisfy: 

3e>0, E|A"| 4+e <oo. 

Assumption A-2. The complex function a : [0, 1] x [0, 1] — ► C is such that \a\ 2 is continuous 
and therefore there exist a non-negative constant cr max such that 

V(x,y)e[0,l] 2 , 0<|a(.T,y)| 2 < ( 7 2 nax <^- (5-7) 

If A„ is a complex deterministic N xn matrix whose non-diagonal entries are zero, assume 
that: 

Assumption A-3. There exists a probability measure H( du, dX) over the set [0, 1] xl with 
compact support TL such that 

1 N 

i— 1 
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In Thcorcm l2 . 31 ( Eg . (|2.6|) and (|2.7|) 'l . one must replace a by its module \a\. The statements 
of Theorem 12.41 and Corollary 12 . 51 are not modified. 
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